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ABSTRACT 

The evolution of quasar clustering is investigated with a new sample of 388 quasars 
with 0.3 < z < 2.2, B < 20.5 and M B < -23, selected over an area of 24.6 deg 2 
in the South Galactic Pole. Assuming a two-point correlation function of the form 
£(r) = (r/r )~ L8 , we detect clustering with tq = 6.2 ± 1.6 h^ 1 comoving Mpc and 
£(r = 15 h~ 1 Mpc) = ^ Jq x 2 ^(x)dx = 0.52 ± 20 at an average redshift of < z = 1.3 >. 
We find a 2a significant increase of the quasar clustering between z = 0.95 and z = 1.8, 
independent of the quasar absolute magnitude and inconsistent with recent evidence 
on the evolution of galaxy clustering. If other quasar samples are added (resulting 
in a total data-set of 737 quasars) the increase of the quasar clustering is still favoured 
although it becomes less significant. With a parameterization of the evolution of the type 
£(r, z) = (r/ro)~ 7 (l + z)^^" /+e ^ we find e ~ —2.5. Evolutionary parameters e > 0.0 are 
excluded at a 0.3% probability level, to be compared with e ~ 0.8 found for galaxies. 
The observed clustering properties appear qualitatively consistent with a scenario of 
Q = 1 CDM in which a) the difference between the quasar and the galaxy clustering 
can be explained as a difference in the effective bias and redshift distributions, and b) 
the quasars, with a lifetime of t ~ 10 8 yr, sparsely sample halos of mass greater than 
M min ~ 10 12 - 10 13 h- 1 Mq. 

We discuss also the possibility that the observed change in the quasar clustering is due 
to an increase in the fraction of early-type galaxies as quasar hosts at high z. 

Subject headings: Quasars: general — Large-scale structure of the universe 
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1. Introduction 

The first detections of the quasar clustering date 
back more than one decade (|Shanks et al. 1983| , 



shaver 1984|) . Up to now, however, more detailed 
studies of the clustering dependence on physical pa- 
rameters like absolute magnitude and redshift was 
hampered by the small number of quasars in statisti- 
cally well-defined s amples. Recent studies of com p lete 
quasars samples ( Andreani and Cristiani 1992 , Mo 
md Fang 1993 , Croom and Shanks 1996 ) found either 



no [ significant difference in the clustering amplitude 
me aaurcd in comoving coordinates at low (a < 1.4) 

and high (z > 1.4) redshift, or marginal detections of 
a decrease of the quasar clustering with redshift. But 
the uncertainties were still quite large. 

The spatial two-point correlation function, £(r), 
expresses the deviation of the observed spatial dis- 
tribution from a random distribution. An estimator 
for the correlation function is £(r) = N°/Np — 1. N° 
is the number of quasars pairs observed at the sepa- 
rations in the range r — Sr/2 to r + dr/2, and N£ is 
the number of pairs counted at the same separation 
in an artificial random catalogue. A common param- 
eterization used to describe the evolving correlation 
function of galaxies is£(r, z) — (r/r )~ 7 (l + z)~( 3+e \ 
where the length is measured in physical (proper) co- 
ordinates and e is an arbitrary fitting parameter. If 
the separation length is measured in comoving coor- 
dinates the above relation becomes 



(1) 



Another estimate of the clustering properties can 
be obtained via the integrated pair counts 3(r) = 
1 + £(r), where £(r) = ^ x 2 ^{x)dx. The advan- 
tage in calculating the integrated two-point correla- 
tion function is that it is stable and does not depend 
on the binning when r is large. When the boundary 

„ N°(<r) 

effect is negligible, H(r) is given by H(r) = N p r ^ <r ^ 
where N°(< r) and N*{< r) are the number of ob- 
served and random simulated pairs respectively, at a 
distance less than r. We assume qo — 0.5, H — 100 
Km s _1 Mpc -1 and comoving distances throughout. 

2. Data 

In order to improve the S /N in the estimate of the 
clustering evolution, we have built a new sample of 
388 quasars down to Bj — 20.5 over a contiguous 



area of 24.6 deg 2 with redshift in the range 0.3 < z < 
2.2. We define quasars as broad emission line AGNs 
having Mb < —23 with qo = 0.5 and Hq = 50 Km 
s _1 Mpc -1 . This is the largest sample of quasars of 
this type yet made. The candidates were observed 
with the Meudon-ESO Fibre Optic System (MEFOS, 



Bellenger et al. 1991) at the ESO 3.6m telescope. 



The sample is centered in the South Galactic Pole 
region where some of the Durham/AAT sample areas 
are included (Boyle et al. 1990), and part of the high 



redshift quasar survey of Warren, Hewett and Osmer 



(1991) was carried out. 206 quasars are newly iden- 



tified. The field is part of the Homogeneous Bright 



Quasar Survey (Cristiani et al. 1995). The quasars 
were selected through the ultraviolet excess (UVx) 
color criterion in the U — Bj versus Bj —R plane. The 
sample is divided in two regions: in the central region 
A (with an area of 9 deg 2 ), including 176 quasars, 
the sample is complete (the surface density is 19.7 
quasars/deg 2 for B < 20.5), while in region B (with 
an area of 16 deg 2 ), including 212 quasars, a fraction 
of the quasar candidates is still not observed, result- 
ing in an incompleteness of about 30%. The sample 
with the full quasar catalogue will be described in a 
forthcoming paper (La Franca et al. , in preparation) . 

3. Methods and Results 

Because of the incompleteness in the region B of 
the quasar sample, the overall nonuniform spatial dis- 
tribution of the quasars complicates the generation of 
the "random" data sets used in the computation of 
the correlation function. In order to guarantee that 
the random data set had exactly the same spatial se- 
lection function as the quasar sample, each of the ran- 
dom data sets was drawn from exactly the same (a, 
S) positions as the overall quasar sample. In each 
random data set, these positions were assigned dif- 
ferent, random, redshifts. We adopted three different 
methods for assigning the random redshifts and im- 
plemented these in three independent codes for com- 
puting £(r). In the three methods the redshifts were 
generated (1) random, from the observed N(z) binned 
in intervals of 0.1 in z\ (2) random, according to the 
distribution expected on the basis of the quasar Lu- 



minosity Function (La Franca and Cristiani 1997); 
(3) scrambling the observed redshifts of the quasars 
in the sample. In the simulation of the quasar red- 
shifts the sample was splitted into two sub-samples: 
a bright (B < 18.7) and a faint (18.7 < B < 20.5) 
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one. This differentiation is necessary as the bright 
sample is complete over the whole (A + B) area (see 
Gristiani ct al. 1995), while the faint sample is the 



result of several MEFOS pointings which are com- 
plete down to B — 20.5 only over the central region 
A, but are only locally complete over the remaining 
region B. Each random catalogue was generated 1000 
times. The results from these three methods agree to 
within the nominal uncertainties obtained from Pois- 
son statistics. 

The data set was divided into several luminosity, 
redshift and spatial sub-samples in order to study the 
autocorrelation function £(r) and the integral auto- 
correlation function £(r) as a function of the comov- 
ing distance, and the correlation length r assuming 
a fixed value of 7 = 1.8. The results are summarized 
in Table 1. In Figure la we show the autocorrela- 
tion function £(r) in the range 0.3 < z < 2.2 for the 
total data set of 388 quasars. One sigma Poisson er- 
rors in the figure are based on the number of observed 
pairs in each bin. The average redshift is 1.34. With 
the slope fixed to 7 = 1.8 it results ro = 6.2 ± 1.6 
h,- 1 Mpc and £(25, 1.34) = 0.21 ± 0.16, in agree- 
ment with estimate of Croom and Shanks (1996) of 
£(25) = 0.16 ± 0.08 (their Table 1). 

In order to examine the evolution in the ampli- 
tude of the correlation function, the sample was split 
into the two redshift ranges 0.3 < z < 1.4, and 
1.4 < z < 2.2, with average redshift of 0.97 and 1.82 
respectively. The resulting £(r) are shown in Figure 
1&. These were fitted by 7 = 1.8 power laws with ro 
as a free parameter. At low redshift ro = 4.2 h^ 1 Mpc 
was found, corresponding to £(15, 0.97) = 0.26 ±0.27; 
while at high redshift ro = 9.1 hr 1 Mpc, which cor- 
responds to £(15, 1.82) = 1.03 ± 0.36, a 1.7a signifi- 
cant discrepancy. The number of observed pairs with 
r < 15 hr 1 Mpc in the range 0.3 < z < 1.4 is 22 
against 17.5 pairs expected. In the redshift range 
1.4 < z < 2.2, 26 pairs are observed against 11.3 
expected. 

In order to check the reliability of this result the 
same kind of analysis was applied in the central com- 
plete region A of our data set, using random samples 
generated by the three methods described above with 
the addition of a random generation of a and 6 coor- 
dinates. In the low redshift subsample, containing 99 
quasars, it turns out £(15,0.95) = 0.12 ± 0.30, while 
in the high redshift subsample, containing 77 quasars 
£(15,1.76) = 1.38 ± 0.48. Also in this case the two 
subsamples show an increase (2.1er significant) of the 
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Fig. 1. — a) The autocorrelation function £(r) for the 
total 388 quasars in the SGP sample with 0.3 < z < 
2.2 as a function of the r comoving distance (q = 0.5, 
H a = 100). The dotted line is the fit for r = 6.2 with 
fixed 7 = 1.8. b) The correlation function £(r) for the 
388 quasars in the SGP sample in two redshift ranges 
0.3 < z < 1.4, and 1.4 < z < 2.2. The dotted line is 
the best low- z fit with ro = 4.2 while the dashed line 
is the best high-z fit with ro = 9.1. c) The correlation 
function £(r) for the complete data-set (737 quasars). 
The dotted line is the best low-z fit with r = 6.1 
while the dashed line is the best high-z fit with ro = 
8.0. 
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Table 1 

Autocorrelation function values with fixed 7 = 1.8 



Area 


Redshift 


ran 


ge 


Magnitude 


range 


Nqso 


<M B 


> 


< z > 




£"(15) 




ro 


A+B 


0.3 < z 


< 1 


4 


M B 


< 


-23 


221 


-25 





0.97 





26 ±0 


27 


4.2+ 2.5 


A+B 


0.3 < z 


< 2 


2 


M b 


< 


-23 


388 


-25 


6 


1.34 





52 ±0 


20 


6.2+ 1.6 


A+B 


1.4 < z 


< 2 


2 


M b 


< 


-23 


167 


-26 


5 


1.82 


1 


03 ±0 


36 


9.1+ 2.0 


A 


0.3 < z 


< 1 


4 


M b 


< 


-23 


99 


-24 


9 


0.95 





12 ±0 


30 


2.8+ 4.9 


A 


0.3 < z 


< 2 


2 


M b 


< 


-23 


176 


-25 


5 


1.30 





74 ±0 


34 


7.1+ 1.9 


A 


1.4 < z 


< 2 


2 


M b 


< 


-23 


77 


-26 


3 


1.76 


1 


38 ±0 


48 


10.8+ 2.6 


A+B 


0.3 < z 


< 1 


4 


-27 < M B 


< 


-25 


106 


-25 


7 


1.11 





09 ±0 


31 


2.3 + 7.5 


A+B 


0.3 < z 


< 2 


2 


-27 < M B 


< 


-25 


231 


-26 





1.47 





82 ±0 


34 


8.0 + 2.2 


A+B 


1.4 < z 


< 2 


2 


-27 < M B 


< 


-25 


125 


-26 


2 


1.77 


1 


10 ±0 


42 


9.5 + 2.6 



clustering amplitude with redshift. 

A possible bias in this type of analysis could be 
introduced by the presence of large group of quasars. 
For this reason we have used the friend of friend tech- 
nique to look for groups containing quasars having at 
least another quasar at a distance less than 15 h~ l 
Mpc. The technique consists of drawing a sphere 
of radius r around each quasar. If there are other 
friends within this sphere they are considered to be- 
long to the same group. The same procedure is ap- 
plied to the friends until no more friends are found 
(see Romberg, Kravtov and Lukash 1996 for details 
about the method). At redshift less than 1.4 there 
are four groups of three members each, while at red- 
shift higher than 1.4 there are four groups of three 
members and one of four members. The four mem- 
ber group is responsible for only 3 of the total 26 
observed pairs at distance less 15 h^ 1 Mpc at high 
redshift, while 11.3 are randomly expected. We can, 
thus, conclude that the result is not contaminated by 
the presence in the sample of a rare large group of 
quasars at z > 1.4. 

We have also investigated possible biases due to 
selection effects in the sample as a function of red- 
shift. Being a flux limited sample, the high redshift 
sample selects preferentially brighter quasars. Natu- 
ral bias predicts that the amplitude of the autocor- 
relation function would increase with galaxy bright- 
ness. The high redshift sample has an average abso- 



lute magnitude of Mb — —26.5, while the low red- 
shift one has Mb — —25.0. In order to disentan- 
gle the magnitude dependence we have computed £ 
for the two subsamples with absolute magnitude in 
the range —27 < Mb < —25. In this case the av- 
erage absolute magnitude for the high redshift sam- 
ple is M B = -26.2 with £(15,1.77) = 1.10 ± 0.42, 
while for the low redshift sample is Mb = —25.7 with 
£(15, 1.11) = 0.09 ± 0.31. Also in this case the two 
subsamples show an increase of the clustering ampli- 
tude with redshift, showing that there is no significant 
dependence of the quasar clustering on the absolute 
magnitude. 

4. Comparison with other samples 



Our clustering result has been compared with the 
cluste ring computed from other complete quasar sur- 



veys (Boyle ct al. 1990, 


La Franca, Cristiani anq 


Barbieri 1992, 


Zitclli ct al. 1992). We excluded the 



fields included in the SGP from the Durham/ A AT 
survey of Boyle et al. (1990). A total of 349 quasars 
have been collected. As already discussed by An- 
dreani and Cristiani (1992), these samples show no 
significant evolution in comoving coordinates. How- 
ever, by adding our new sample to these data (a total 
data-set of 737 QSOs) the statistical significancy on 
the increase of the quasar clustering evolution with 
redshift is reduced. The resulting £(r) are shown 
in Figure lc. At low redshift we find £(15,0.97) = 
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Fig. 2.— The amplitude of the £(15 h^ 1 Mpc) as 
a function of redshift. The points are the low and 
high redshift SGP subsamplcs (filled circles); our SGP 
sample plus the high redshift complete quasar sam- 
ples from Boyle et al. (1990), La Franca, Cristiani 
and Barbieri (1992), and Zitelli et al (1992): complete 
sample (open circle), and divided in the two redshift 
subsamples (open triangles); the low redshift AGN 
(Seyfertl and Seyfert2) amplitude from Boyle and Mo 
(1993) and Georgantopoulos and Shanks (1994) (filled 
triangle); the high redshift sample from Kundic 1997 
(open square). The dotted line is the e = —2.5 clus- 
tering evolution, while the dashed lines are the 10 12 
and 10 13 Mq h~ x minimum halo masses clustering 
evolution according to the transient model of Matar- 
rese et al. (1997). 



0.50 ±0.24 (r = 6.1 hr 1 Mpc), while at high redshift 
£(15,1.85) = 0.80 ±0.28 (r = 8.0 h' 1 Mpc). The 
total number of observed pairs with r < 15 h^ 1 Mpc 
in the range 0.3 < z < 1.4 is 39 against 26 pairs ex- 
pected. In the redshift range 1.4 < z < 2.2 a total of 
42 pairs are observed against 23.3 expected. In the 
complete redshift range 0.3 < z < 2.2 we find an av- 
erage value of £(15, 1.37) = 0.64 ± 0.18 corresponding 
to r = 7.0 h~ x Mpc. 

In Figure 2 the evolution of £(15) as a function of 
redshift is shown. We compare our result with previ- 
ous analyses of the quasar clustering at low and high 
redshift. Boyle and Mo (1993) obtained clustering 
statistics for low redshift (z < 0.2) quasars in the Ein- 
stein Extended Medium Sensitivity Survey (EMSS). 
Below 10 h^ 1 Mpc four quasar pairs are observed 
compa red with the prediction of 2.37 pairs ( £(10) 
i i. 7 : Georgantopoulos and Shanks (1994)| carried 



out a similar analysis on the low redshift (z < 0.1) 
Seyferts 1 and 2 galaxies in the IRAS Point Source 
Catalogue finding £(20) = 0.52±0.13. For the Seyfert 
1 only sample it turns out £(20) = — 0.10±0.27. In or- 
der to compare these results with our data, we added 
together the EMSS and IRAS samples and extrapo- 
lated the numbers of pairs below 10 and 20 hr 1 Mpc 
respectively to a distance of 15 h^ 1 Mpc, assuming 
a power law slope 7 = 1.8. This point is shown in 
Figure 2 with £(15, 0.05) = 0.24 ± 0.25. 

We have fitted the two points at low and high red- 
shift (obtained from our sample plus the quasar sam- 
ples from Boyle et al. (1990), La Franca, Cristiani 
and Barbieri (1992), and Zitelli et al. (1992)) and 
the point from the EMSS and IRAS sample with the 
£ dependence on redshift expressed in eq. 1, assum- 
ing 7 = 1.8. We find e = -2.5 ± 1.0 with £(15, 0) = 
0.22 ±0.20, which corresponds to r (z = 0) = 3.9 hr Y 
Mpc or r (z = 1.5) = 7.3 ft -1 . From x 2 statistics the 
two values of £(15) in the redshift range 0.3 < z < 2.2 
are indistinguishable from a constant comoving clus- 
tering evolution. However, statistics are sufficient to 
constrain the quasar clustering evolution. All the evo- 
lution models normalized to the £(15,0.93) = 0.51 
average quasar and local Seyfert clustering and with 
e > 0.0 are rejected at a probability level less than 
0.003, while the models with a normalization to the 
£(15, 1.37) = 0.64 average quasar clustering and with 
e > —0.5 are rejected at a probability level less than 
0.003. As discussed in the next section, this result 
should be compared with e ~ 0.8 found for galaxies. 

A recent analysis of the quasar correlation func- 
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tion in the Palomar Transit Grism Survey has found 
that £0 > 2)/£{z < 2) = 1.8±?;1 (Kundic 1997), 
with an average redshift of 1.3 for the low-redshift 
subsample and 3.1 for the high redshift one. The 
corresponding high redshift point (£(15,3.1) = 1.2) 
is plotted in Figure 2 and was obtained by multi- 
plying our low-redshift £(r = 15) measurement at 
z = 1.37 by i{z > 2)/£(z < 2) given by |Kundic 
(1937). Stephens et al. (1997) from an indcpcn- 



contribute already at ~ 28% of the total XRB (see 



McHardy et al. 1997 their Figure 7) and this corre- 



dent analysis of the same data-set find an even higher 
value with r$(z = 3.3) = 17.5 h~ x corresponding to 
£(15,3.3) = 3.3. These analyses, although uncertain- 
ties on the clustering of this high redshift sample are 
very large, strengthen the evidence for an increase of 
the quasars clustering amplitude with redshift and are 
in agreement with the e = —2.5 evolution law within 
a lcr confidence level. 

5. Limits imposed by the XRB ACF 

Studies of the X-ray background (XRB) residual 
fluctuations and their autocorrelation function (ACF) 
have been used to constrain the counts, evolution 
and clustering properties of the extragalactic X-ray 
sources, in particular quasars. Higher values of Tq 
and lower values of e determine a larger contribution 
of quasars to the soft XRB fluctuations, with the pos- 
sible result of exceeding the limits on the XRB ACF 



and spectral shape. The general consensus (e.g. Car 
:era and Barcons 1992; Danese et al. 1993, Soltan and 



Hasinger 1994 ) is that a population of sources, clus- 
tered like quasars, cannot produce more than 30% - 
50% of the still unresolved XRB. Soltan and Hasinger 
1994 show (their Figure 8) that if quasars contribute 
100% of the total XRB they would require a cor- 
relation length of 1.5-2 hr 1 Mpc at z = 1.5 with 
— 1.2 < e < 2.4. If instead a quasar contribution 
of 35% is assumed, a correlation length of 4.5-6 hr 1 
Mpc is allowed. 



A recent analysis of deep ROSAT fields (McHardy 



t al. 1997) has shown that quasars (i.e. broad emis- 



sion line AGNs) unlikely contribute to more than 40 
per cent (but at least 30%) of the total XRB at 1 
KeV. In the framework of this result, our outcomes 
agree with the limits estimated by Danese et al. 1992| . 
Indeed, from Danese et al. (1993, their Figure 3), it 
is possible to see that our estimate of tq(z — 0) = 3.9 
h^ 1 Mpc and e ~ —2.5 allows a contribution of more 
than 15% for fluxes fainter than 5 x 10~ 15 erg cm~ 2 
s _1 in the 0.9-2.4 KeV band. At these fluxes quasars 



sponds to an allowed contribution of quasars to the 
total XRB of more than 45%, in agreement with the 
limits imposed by McHardy et al. (1997). 

6. Discussion 

An increase of the clustering amplitude with red- 
shift is a powerful diagnostic for models of quasar 
formation and evolution. Schematic patterns of the 
type "continuous-activity", in which each AGN un- 
dergoes a long, continuous dimming after a short for- 
mation phase at z > 2, and "recurrent-activity", in 
which the potential- AGN population shows intermit- 
tent episodes of brightening lasting a relatively short 
time, cannot be easily reconciled with the present ob- 
servations. If the same population of mass condensa- 
tions, formed at a given redshift, is permanently or 
recurrently shining as quasars, then in any hierarchi- 
cal clustering scenario an increase of the clustering 
with decreasing redshift would be expected. On the 
other hand, the quasar phenomenon may represent a 
single, short event (t ~ 10 8 yr) in the host galaxy life- 
time. If we associate it with a characteristic mass at 
all epochs, i.e. a sparse sampling of the distribution 
of halos with mass greater than M m i n , quasars would 
correspond to rarer and rarer over-densities with in- 
creasing redshift and their clustering amplitude is ex- 
pected to grow with z. This scenario is described as 



the transient model by Matarrese et al. 1997 . 

In a more quantitative way we can follow Matar- 
rese et al. (1997) in estimating the observed correla- 
tion function in a given redshift interval in the frame- 
work of an = 1 CDM model. This is carried out 
by suitably weighting the mass autocorrelation func- 
tion £(r, z) with the number of objects as a function 
of redshift and convolving with the effective bias fac- 
tor b e ff. The effective bias b e ff(z) is the weighted 
mean of the bias b(M, z) with the mass distribution 
(e.g. Press-Schecter) at a given redshift. Figure 2 
shows that for a minimum halo mass in the range 
M mm = 10 12 - 1O 13 /i _1 M , the QSO two-point cor- 
relation function is reproduced both in absolute value 
and redshift evolution. It is interesting to note that 



Cavaliere et al. 1997, on the basis of the shape and 



evolution of the QSO Luminosity Function, suggest 
that the quasar phenomenon at z < 3 is connected 
with interactions in relatively small galaxy groups, 
whose typical mass, 5- 10 12 Mq, is remarkably similar 



G 



to our estimate, based on the clustering properties. 

What is the relation between the possible observed 
increase of the quasar clustering with redshift, the re- 
cent evidence of a decrease of the galaxy clustering 



with redshift with e ~ 


> 0.8 (e.g. 


Le Fevre et al. 1996), 


Carlberg et al. 1997 




Villumscn, Freudling and da 



Coiita 1997) and the theories of formation and evolu- 
tion of cosmic structure? 

We can schematically envisage that either galax- 
ies form via the merging of lower-mass halos or they 
form at some characteristic redshift (the merging and 
object- conserving models of Matarrese et al. 1997, 
respectively). In this framework the non-linear pro- 
cesses leading to the galaxy formation give origin to a 
bias parameter with a distinctive evolution. The con- 
volution of the bias b(M, z) with the mass autocor- 
relation function £(r, z) and the galaxy redshift dis- 
tribution can produce the observed decrease of the 
galaxy clustering with increasing redshift. 

It must be pointed out, however, that different 
population of galaxies may show different clustering 
properties and there are indications that for high-z 
star- forming galaxies ( Steidel et al. 1997 ) the clus- 
tering amplitude do es not decrease so qu ickly as for 
faint low-z galaxies ( Le Fevre et al. 1996|) . 

Another (not necessarily alternative) scenario, which 
can explain the QSO clustering behaviour, is based 
on a change of the morphological mix of the quasar 
host galaxies with redshift. Recently Franceschini and 



Gr atton (1997) have proposed, on the basis of the 
metallicities and abundance ratios in absorption-line 
systems associated with quasars, that the host galax- 
ies at high redshift are early- type bulge dominated, 
formed quickly through a huge episode of star forma- 
tion. This picture is reinforced by the observed in- 
crease with redshift of the dust mass associated with 



quasars (Andreani, Franceschini and Granato 1997). 
On the other hand, low-z quasars may be related to 
the accretion onto black holes in the nuclei of late- 
type galaxies, where gas is still available. At low red- 
shift early-type galaxies are observed to cluster more 
strongly than late type galaxies (Loveday et al. 199£ ) 
by a factor 3 — 5. If this behaviour persists at all red- 
shifts, an increase of the fraction of elliptical galax- 
ies as quasar hosts could produce the increase of the 
quasar correlation at high rcdshifts. 
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